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ABSTRACT 

The  probability  function  expressing  the  fractional  kill  of  a 
circular  target  engaged  by  a  salvo  of  area  kill  weapons  can  be  expressed 
in  terms  of  a  finite  series.  When  the  salvo  contains  many  rounds, 
the  nature  of  the  series  makes  it  useless  in  computing  the  kill 
probability.  In  cases  where  the  distribution  of  hits  is  bivariate 
normal,  with  equal  or  near  equal  standard  deviations,  this  difficulty 
may  be  overcome  by  expressing  the  series  in  terms  of  a  definite  integral, 
readily  solved  by  numerical  quadrature.  Additionally,  the  definite 
Integral  is  expressed  in  terms  of  a  different  finite  series  having 
properties  making  it  easily  sunsnable  on  a  digital  computer. 
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I.  INTRODUCTION 


An  important  problem  in  weapon  systems  analysis  is  the  computation 
of  the  expected  fractional  kill  of  a  circular  target  when  engaged  by  a 
salvo  of  area  kill  weapons.  The  solution  of  this  problem,  f^,  usually 
is  in  the  form  of  a  finite  series 

N 

fN  =  C  £  (-l)i+1(”)[xVi]  g(i),  (1) 

i=l 

where  C  and  x  are  constants,  0  <  x  <  1,  N  is  the  number  of  rounds  in 

the  salvo,  and  g(i)  is  a  function  dependent  on  the  summation  index  i. 

1* 

This  series  is  developed  by  Groves  and  by  Dr.  F.  E.  Grubbs  in  an 

unpublished  work.  This  sum  can  usually  be  obtained  quite  easily  on  a 

digital  computer;  however,  for  large  N  the  nature  of  the  series  makes 

it  useless  in  computing  f  .  In  cases  where  the  distribution  of  hits 

is  bivariate  normal  with  equal  or  near  equal  standard  deviations,  g(i) 

can  be  expressed  in  the  form 

g(i)  =  (l-e“U),  (2) 

where  X  is  a  constant.  For  this  form  of  g(i),  an  alternate  form  for 

f„  can  be  obtained  which  is  more  amenable  to  computation.  The 

2  3 

procedure  developed  by  Wheslon  leads  to  a  definite  integral  which 
can  be  expressed  in  terms  of  a  finite  series.  This  series,  in  contrast 
to  the  original,  can  easily  be  summed  for  large  N  on  a  digital  computer. 

The  procedure  consists  of  equating  pert  of  the  general  term  of  the 
series,  H(i),  with  the  Laplace  transform  of  a  function  h(y).  Let 

H(i)  -  (l-e‘Xi)/i  -  J  e’^X"1  {(l-e'Xi)/i)  dy,  (3) 

o 


Superscript  numbers  denote  references  which  may  be  found  on  page  11. 

The  nature  of  this  computational  difficulty  and  the  presentation  of 
a  method  to  overcome  it  for  a  more  general  problem  than  the  one 
solved  here  are  contained  in  a  forthcoming  report  by  the  author  and 
Lynn  S.  Mohler. 
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where  (l-e"’^i)/i}  is  the  inverse  Laplace  transform  of  (l-e  ^J/i, 
and  i  is  treated  as  the  transform  variable.  Now 

£_1t<l-e'X1>/i)  -  Cl-Sx(yn,  (U) 

where  S^(y)  is  the  Heaviside  step  function,  i.e., 

* 

0,  0  <  y  <  X, 

S.(y)  -  (5) 

1,  y  >  X. 

If  we  let  FN  =  fN/C,  may  be  expressed  in  the  form 

fn=  I  <-1)1+1(i)xl  Jo  «'ly  Ci-Sxtyn  ay 

=  J  tl-\(y)1  Y,  <-1)i+1(  J)  *ie”ly  ay.  (6) 

°  i=*l 

Letting  w  =  xe  y,  we  note  that 

i-l  1=1 

«  [l-(l-w)N].  (T) 

Hence 

FN-f  [l-Sx(y)][l-(l-xe"y)N]  dy 

0 

.  J  U-d-xe^)11]  dy  ,  (8) 
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(9) 


fn  =  x  "  JX  (1-xe’y)N  <*y. 

o 

Due  to  the  form  of  the  result  in  Equation  (?),  it  becomes  apparent  that 
proper  substitutions  in  the  original  multiple  integral^  and  the  use 
of  appropriate  transformations  of  variables  would  lead  to  the  result 
of  Equation  (9)*  Nevertheless  the  approach  used  seems  instructive  in 
itself. 


II.  AN  ALTERNATE  SERIES 
If  in  Equation  (9)  we  let 

u  =  (l-xe"y),  (10) 

we  obtain 


where  u^ 


(l-x)  and  u^ 


(1-xe  X). 


Since 


(l-u)  =  £  u1,  0  <  u  <  1, 


i=0 


F„  can  be  written  as 
N 


F„  -  X 


I  r « 

i-0  ul 


i+N 


du 


X  -  l  (uX-u^)/i  . 


1  “N+l 


(11) 


(12) 


1 

i 

1 


Q 


Since 


N 

=  £  [(l-xe-y  -  (l-x)l]/i.  (1U) 

1=1 

The  terms  of  the  series  in  Equation  (iU)  are  well  behaved,  i.e.,  they 
decrease  monotonically  with  the  summation  index  i  and  are  therefore 
readily  sumraable. 

This  series  should  be  useful  in  computing  kill  probability  for 
large  N  where  g(i)  can  be  expressed  by  Equation  (2).  Additionally, 
it  has  proved  useful  in  checking  the  accuracy  of  more  general  methods 
used  to  solve  the  problem  in  cases  where  g(i)  cannot  be  expressed  by 
Equation  (2). 
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